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PREFACE 


Compact  microwave  antennas  formed  from  resonant  microstrip  elements 
have  recently  attracted  considerable  attention  because  of  their 
simplicity  and  conformability  to  shaped  surfaces. 

The  Department  of  Transportation  has  been  especially  interested 
in  this  type  of  antenna  for  possible  application  with  the  Global 
Positioning  System  (GPS)  and  in  aircraft  satellite  communications 
systems. 

The  work  reported  in  both  Volumes  I  and  II  of  this  document  is 
aimed  at  providing  a  general  theoretical  framework  useful  for  analyz¬ 
ing  and/or  synthesizing  microstrip  antennas  and  arrays.  The  work  was 
completed  and  this  report  prepared  by  F.  R.  Morgenthaler  under  the 
direction  of  TSC  Technical  Monitor  Leslie  Klein. 
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EXECUTIVE  SUMMARY 


Microwave  antennas  formed  from  resonant  microstrip  elements  have  recently 
>  attracted  considerable  attention  because  of  their  conformability  to  shaped  sur¬ 
faces,  and  their  compactness  and  simplicity. 

Although  calculations  based  upon  relatively  simple  theoretical  models 
for  the  microstrip  resonators  appear  to  be  in  reasonable  agreement  with  the 
experimental  results  available  for  comparison,  it  is  of  interest  and  benefit 
to  critique  the  approximations  that  ''ive  been  employed.  This  is  so,  because 
high  performance  from  an  array  of  elements  often  can  be  difficult  to  achieve 
unless  higher-order  terms  pertinent  to  a  single  element,  and  mutual  couplings 
among  elements,  are  properly  understood. 

With  this  goal  in  mind,  we  first  present  a  general  approach  that,  in 
principle,  would  give  a  complete  description  of  the  radiation  fields  from  any 
•  microstrip  antenna  consisting  of  one  or  more  perfectly  conducting  metal  patches, 
deposited  upon  the  surface  of  a  dielectric  slab  backed  by  a  perfectly  conducting  metal 

f 

sheet.  Certain  of  the  patches  are  presumed  resonant  at  microwave  frequencies 
of  interest. 

In  particular,  we  review  and  then  employ  an  important  equivalence 
theorem  that  allows  the  radiation  1 d s  to  be  calculated  from  the  components 

of  electric  and  magnetic  fields  that  are  tangential  to  any  surface  enclosing 
all  of  the  radiating  currents  and  charges. 

The  approach  used  is  first  to  find  the  modal  structure  within  this 
surface  (as  effected  by  assumed  fringing  and  radiation  of  the  outer  fields) 
and  second  to  calculate  these  effects  in  a  self-consistent  manner. 

Specifically,  we  develop  formulas  to  estimate  the  strength  of  the 
electric  charges  and  electric  currents  that  are  present  on  the  top  of  the 
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structure  consisting  of  coupled  X/4-resonant  microstrip  patches,  the  dielec¬ 
tric  layer,  and  the  ground  plane,  and  the  components  at  the  radiation  field 
arising  from  the  currents. 

In  order  to  accomodate  dual -frequency  operation.  Ball  Aerospace  Systems 
Division  has  developed  a  stacked  element  concept  in  which  the  smaller  high- 
frequency  patch  is  placed  above  the  larger  low-frequency  one. 

Because  this  is  an  important  concept,  we  develop  a  transmission  line 
model  of  this  element  with  distributed  coupling  and  solve  for  the 
excitations  in  both  patches  when  the  drive  frequency  corresponds  to 
either  the  low  or  high  resonance  condition.  We  consider  the  cases  of 
substrates  both  with  dissimilar  and  with  identical  dielectric  properties. 

Throughout,  emphasis  is  placed  on  developing  equations  describing  both 
the  near  fields  and  the  radiation  fields;  from  the  latter  it  is  a  simple  mat¬ 
ter  to  calculate  radiation  conductance,  efficiency  and  pattern. 

As  an  aid  in  these  endeavors, we  develop, in  the  Appendix  A,  a  novel  approx¬ 
imation  that  helps  provide  physical  insight  with  respect  to  field  patterns, 
coupling  between  patches  and  the  like.  For  rectangular  microstrip  patches, 
the  approximation  allows  simple  analytic  expressions  to  be  obtained  for  the 
electric  and  magnetic  field  coupling  coefficients  between  patches. 

As  an  important  by-product,  several  new  approximate  formulas  are  ob¬ 
tained  that  very  accurately  predict  the  characteristic  impedance  and  phase 
velocity  of  microstrip  transmission  lines  of  arbitrary  width  when  the  dielec¬ 
tric  constant  of  the  substrate  is  also  arbitrary.  A  comparison  with  other 
formulas  established  in  the  literature  reveals  that  the  relative  accuracy  is 
within  1-2  percent. 

In  an  effort  to  highlight  the  fundamental  formulation  and  facilitate 
use  of  the  key  results,  especially  important  equations  have  been  placed  within 
a  rectangular  box  outline. 
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1.  INTRODUCTION 


Microwave  antennas  formed  from  resonant  microstrip  elements  have  recently 
attracted  considerable  attention  because  of  their  conformabil ity  to  shaped  sur¬ 
faces,  and  their  compactness  and  simplicity. 

1  2  3 

The  early  work  of  Deschamps  was  extended  by  Munson,  Itoh  and  Mittra, 

4  5  6 

Sanford,  Weinschel,  and  more  recently  elaborated  upon  by  Sanford  and  Klein, 

7  8  9  10  11 

Howell,  Agrawal  et  al Derneryd,  Lo  et  al . ,  James  et  al.  and  others. 

Although  calculations  based  upon  relatively  simple  theoretical  models 
for  the  microstrip  resonators  appear  to  be  in  reasonable  agreement  with  the  ex¬ 
perimental  results  available  for  comparison,  it  is  of  interest  and  benefit  to 
critique  the  approximations  that  have  been  employed. 

This  is  so  because  high  performance  from  an  array  of  elements  often 
can  be  difficult  to  achieve  unless  higher-order  terms  pertinent  to  a  single 
element  and  mutual  couplings  among  elements  are  properly  understood. 

With  this  goal  in  mind,  we  first  review  a  general  approach  that,  in 
principle,  would  give  an  exact  description  of  the  radiation  fields  from  any 
microstrip  antenna  array  consisting  of  one  or  more  perfectly  conducting  metal 
patches,  such  as  that  shown  in  Fig.  1,  deposited  upon  the  surface  of  a  dielec¬ 
tric  slab  of  finite  extent  and  of  thickness  h,  and  backed  by  a  perfectly  conducting 
metal  sheet.  The  sandwich  structure  is  surrounded  by  free  space. 

The  basic  characteristics  of  microstrip  lines  are  developed  in  Ap¬ 
pendix  A  by  novel  approximations  that  help  provide  physical  insight,  with 
respect  to  field  patterns,  coupling  between  lines  and  the  like.  For  the 
present  discussion,  such  details  can  be  postponed. 
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GURE  1.  BASIC  GEOMETRY  OF  MICROSTRIF  PATCH 


2 


Although  the  actual  sources  of  E  and  H  consist  of  free  electric 

currents  flowing  on  the  conducting  surfaces,  free  charge  resident  upon  them, 

and  polarization  current  and  charge  on  and  within  the  dielectric  slab,  it  is 

permissible  to  utilize  fictitious  magnetic  current  and  charge  densities 

(J*,  p*)  (in  addition  to  electric  current  and  charge  densities  (J,  p))  to 

equivalently  represent  the  sources  at  points  external  to  them. 

12 

Following  Chu  ,  we  express  Maxwell's  Eqs.  in  MKS  units  as 


v  *  s .  £o  14  .  j 

(1) 

V  x  E  +  Uq  _  =  -J 

(2) 

£0  V-  E  =  p 

(3) 

U0  V-  H  =  p*. 

(4) 

Since  these  equations  are  linear,  it  is  convenient  to  resolve  E  and  H  into 
components  E^ ,  ,  due  to  J,  p  alone  and  E2,  H2  due  to  J*,  p*  alone. 

Then 


Vi0H1  =  V  x  A 

(5a) 

Ei  ■  -It-7* 

(5b) 

where 

A(fp,t)  =  - 

|  v!iv t_  dv 

1  rQp  Q 

V 

(6a) 

<J>(rp,t)  =  - 

|  pira  — j* 1  dv 

1  <"<■„  ■‘op  q 

(6b) 

c  =  ^  is  the  velocity  of  light  in  free  space,  rn  is  the  distance 

*Vo  _  _gp 

between  the  source  point  Q  and  field  point  P,  and  r^  and  rp  the  respective  vector 

distances  as  measured  from  a  common  origin. 
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Analogously, 

where 

and 


At  a  surface  of  electromagnetic  discontinuity  at  least  one  of  the  volume 
densities  3,  p,  3*,  p*  is  infinite  and  the  respective  boundary  conditions  as¬ 
sociated  with  Eqs.  (1)  -  (4)  are 


nx(H^  -  H^)  =  K 

(9a) 

nx(E(])  -  E(2))  =  -R* 

(9b) 

Vo(l(ir  E(2)’  =  °s 

(10a) 

V  %,  -  V)>  ■  °S 

(10b) 

where  the  direction  of  the  normal  n  with  respect  to  the  surface  and  the  de¬ 
finition  of  sides  (1)  and  (2)  are  given  in  Fig.  (2).  Note  that  H^. j  should  not 
be  confused  with  ;  the  number  in  parenthesis  refers  to  the  side  of  the  sur¬ 
face.  Thus  refers  to  the  field  of  Eq.  (5a)  on  side  (2)  of  the  surface 

of  discontinuity. 

Here  K  and  K*  are,  respectively,  electric  and  "magnetic"  surface  cur¬ 
rent  densities  and  o$  and  o*  are,  respectively,  electric  and  "magnetic"  surface 
charge  densities. 


Conventions  pertaining  at  a  surface  of  discontinuity.  Th< 
unit  normal  vector  n  is  directed  from  side  (2)  toward 
side  (1). 

FIGURE  2.  SURFACE  OF  ELECTROMAGNETIC  DISCONTINUITY 


In  a;  the  volume  sources  3,  p,  3*,  p*  within  the  closed 

surface  S  create  E,  H  outside.  In  b,  appropriate  sources  K,  o 

-  *  o  * 

K  ,  s  located  on  the  surface  create  the  identical  fields  out 
side,  but  0-0  and  H  -  0  within  the  surface. 


FIGURE  3.  EQUIVALENCE  PRINCIPLE  FOR  ELECTROMAGNETIC  FIELDS 


stabllsh  a  very  W*** 

These  results  may  be  used  to  es  ^  closed  surface  s 

.....  — 

rr.T— -  ■  r~ 

. .  -  — - — * -  - 

„„  the  surface.  Ihe 

tlcal  to  those  1.  ««•  ^ 


-  n  h  ,  .  0  and  the 

a,  vector  poihts  outward,  t(2)-  ■  (2) 

Since  the  normal  vecto 

required  surface  sources  are  _  ( 


(Via) 


K  =  nxH 

_Y*  =  n  x  E 
0S  *  Eo"'  E 


*  =  y  n  •  H. 

$  purees  entemal .  it 

•ded  only  that  there  are  no 

As  \s  ^-Pr°;;b)  ^  (8a,b)  can  he  replaced  hy 

therefore  /  ....  -  rnn/0  V ^  03a) 


(11b) 

(12a) 

(12b) 


Mfpnt)  -  + 


^(Fp.t)  =  ’ 


A*(rp.t) 


(13a) 


S  l 

=  +  $>  Anr- 


(13b) 

(14a) 
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and 


(14b) 


i|>*(rp,t)  =  - 


H(?Q>t:'rQp/c),ndaQ 


4nr 


QP 


The  element  dag  is  the  scalar  differential  area  at  any  point  Q  that  lies  on 
the  surface  S. 

Notice  that  for  complex  harmonic  time  dependence,  employing  retarded 

t1me  exp  ju>(t-(-3£  ))=  exp(-jkrgp)  exp(joit) 

where  k  =  w/c.  In  this  case,  the  complex  vector  potentials  for  sinusoidal 
;ady  state  frequency  m  are 


It  is  not  necessary  to  calculate  the  scalar  potentials  from  the  normal 
components  of  the  fields,  because,  in  free  space,  the  total  fields  are  expres- 
sible.js 

E  =  -j  c/kVx  (V  x  A)  +  1/eqVx  A*  (lta) 

H  =  j  c/k V x  (V  x  A*)  +  1/p0Vx  (lf>b) 

This  is  a  consequence  of  their  having  to  satisfy  the  complex  form  of  the  free 
space  Maxwell  Equations. 


,  •  ...  :>;> 
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SUMMARY  REMARKS 


Knowledge  of  the  E  and  H  fields  tangential  to  any  surface,  S,  enclosing 
all  currents  and  charges  allows  A  and  A*  to  be  calculated;  these  in  turn  allow 
the  fields  to  be  obtained  at  any  point  external  to  S.  When  that  point  is  far 
from  the  antenna,  even  fairly  crude  approximations  of  the  fields  on  S  will 
allow  the  radiation  fields,  radiation  conductance  and  pattern  to  be  calculated 
in  a  straight-forward  manner. 

This  information,  when  fed  back,  allows  the  near  fields  to  be  calculated 
with  improved  accuracy. 

We  turn  next  to  the  important  case  of  the  rectangular  microstrip  patch 
radiator  -  essentially  a  X/4  strip  transmission  line  radiating  from  the  open 
circuit  end  of  the  resonator.  [We  employ  reasonable  approximations  to  derive 
A  and  A*].  The  element  is  driven  at  an  appropriate  distance  from  the  short 
circuit  end  where  the  impedance  allows  a  favorable  match  to  the  input  micro¬ 
strip  transmission  line. 
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2.  RECTANGULAR  MICROSTRIP  PATCH  RADIATOR 

Consider  a  single  rectangular  patch  of  width  w  and  length  L  shorted  to 
the  ground  plane  at  z  =  0  as  shown  in  Fig.  4a,  b.  Assume  first  that  the  dielec¬ 
tric  constant  of  the  supporting  slab  is  unity. 

If  the  resonator  is  excited  at  the  frequency  for  which  L  is  approximately 
A/4,  the  mode  pattern  will  resemble  that  shown  in  the  figure.  In  particular,  if 
h/A  <<  1,  the  fields  between  the  plates  will  be  quasi-TEM  and  the  fringing 
fields  at  the  edges  will  be  comparatively  weak. 

The  equivalent  circuit  of  the  transmission  line  will  then  be  as  shown 
in  Fig.  5  where  Gq  and  Bq  are  in  general  frequency  dependent  when  c‘>  1 
(because  then  exact  TEM  modes  are  not  possible).  For  the  present  c'  =  1,  there¬ 
fore  Bq  =  0  and  the  frequency  independent  GQ  includes  the  effects  of  lateral 
fringing  (Fig.  4b).  The  lumped  capacitance  models  the  fringing  at  the  open 
circuit  (shown  in  Fig.  4a).  The  conductance  Gr  represents  the  power  radiated 

1  ^ 

from  the  structure  or  else  lost  due  to  conversion  to  surface  wave  excitations. 

Since  normally  the  resonator  is  high  Q,  may  be  neglected  when  solving 
for  the  resonant  frequencies. 

Transverse  resonance  (or  what  is  the  same .continuity  of  voltage  and 
current)  requires 

For  wCf/Yo  <<  1,  the  principal  mode  satisfies  kL=  n/2  +  A 

ioCf 

tan  A  =  A- - y~  •  (18) 

o 

off 

Defining  kL  =  tt/2  and  using 


(19) 


1 


E  and  ft  field  patterns  associated  with  a  resonant  rectangular 
microstrip  patch  of  width  w  and  length  L  *  A/4.  An  rf  short 
exists  at  z  *  0;  the  three  other  sides  are  open. 


FIGURE  4.  E  AND  ft  FIELD  PATTERNS  ASSOCIATED  WITH  A  RESONANT 
RECTANGULAR  MICROSTRIP  PATCH  OF  WIDTH  w  AND  LENGTH 
L  -  A/4 
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cr  _ 

< 

—  < 

.  ""  -  -d 

An  equivalent  transmission  line  model  corresponding  to 
wave  propagation  directed  along  the  z -coordinate  of  Fig.  4 


FIGURE  5.  AN  EQUIVALENT  TRANSMISSION  LINE  MODEL 


gives 


where  C'  is  the  capacitance  per  unit  length  of  the  microstrip  line. 


RADIATION  FIELDS 


When  the  plate  spacing  h  is  very  small,  the  fringing  fields  outside 
of  the  plate  may  often  be  neglected.  In  this  case  then  the  surface  S  of 
Eqs.  (15)  corresponds  to  that  of  the  rectangular  box  with  volume  L  x  W  x  H. 
The  vector  potential  A  = 0  and  A*  contains  contributions  from  the  three  open 
sides. 

For  z  =  -L  | x |< w/2  0  <  y  <  h. 


EMyEo 


(22a) 


and  applying  Eg.  (15b)  gives 


.  £„Ewhe-Jkr  slrfe 
_  .j  o  o _  Zr 

front  face”  x  4nr  ,kwxx 

_  '~2r' 


/JTy 

(7-nr1) 

x  j  y7 


-jkLf 


(22b) 


For  the  sides  x  =  ±W/2,  0  <  y<_  h,  -L  <  z  <  0 


E -  i  E  sinKz 

y  0 


and  again  applying  Eq.  (15b)  gives  r 

£.E.Lhe"^r  / 

—  si 


TXT 


A* 

Hsides 


00 
2ttt” 


kwx  /,^-A  T  ze  "  r  *  J? 
tnfe  (-  li,  )  — - r2 -  (23b) 

(J)2  -  #,  )2  • 


(23a) 


Because  kh«l,  the  factors  involving  y  are  approximately  unity. 

The  principal  effect  of  the  fringing  is  not  to  alter  A*  but  rather 
to  create  a  vector  potential  A  that  is  also  z-directed. 
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The  vector  potential  A  arising  from  the  surface  currents  due  to  the 
fringing  field  can  be  estimated  by  integrating  the  tangential  H-field  along 
the  path  shown  in  Fig.  6  for  a  strip  of  width  dz.  For  a  point  P  in  the 
far  field  region,  Eq.  (15a)  has  only  a  single  component  that  may  be  approxi¬ 
mated  as 


%  e-Jkr 

4n  r 


(  0 


df.dz’H  U  2V^1n0<cos«x-  +  sin^y')  df;dz'  +  cosez' ] 


i  z'=-L 


(24) 


where  ^(either  x'  or  y * )  is  the  surface  coordinate  along  the  path. 
Assuming  a  TEM  field  of  the  form 


Htan(£»*')  c  -Hr(S)  sin(^j-  z ' ) 


K 


and  both  kh  <<1  and  kS  «1,  Eq.  (24)  can  be  further  approximated  as 

-w/2  -  A 

F  Jltr?F«»  1  Hx(x)eJks,oecos«‘x 


w/2  -  A 


Hx(x)ejksinecos*x  dx  + 


Hx(x)ejlcs1n0co^x  dx 


-w/2  +  A 

w/2  +  A 

b 

1 

' 

a 

+  -jk  w/2  sinOcos«|) 

HtU)d£  ♦ 

Ht(c)dF, 

J 

(25) 


where 


F(0) 


- 5 - j  +  - 

(kLcosOT  -  (n/2T  kLcosO  +  n/2 


-j(kLcos0  ♦  n/2)  -j(kLcos0  -  n/2) 

“  C 


kLcose  -  n/2 


(26) 


For  h/w«l,  the  values  of  Hx(x)  for  |x|<w/2-A  and  |x|>w/2  +  A  can  be  ap¬ 
proximated  with  the  help  of  a  vector  potential  Az  that  has  the  same  form  as 
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y 


The  contour  of  integration  associated  with  the 
C  coordinate  of  Eq.  (24). 

FIGURE  6.  THE  CONTOUR  OF  INTEGRATION 
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Eqs.  (A-l )  and  (A-2) of  the  Appendix,  (except  with  e'  =  1). 


The  result  is 


-H  hw 
o 


H  ( x)  ~  .  hw  +  1TC(w/2)^  -xd] 

H  hw 
o 


hw  +  Ti[x  -  (w/2)‘ 


| x  |<w/2  -  A 


j  x  |>  w/2  +  A 


where  Hq  is  the  amplitude  in  the  uniform  field  region  (inside  of  the  parallel 
plates).  These  expressions  allow  the  first  three  integrals  on  the  right  hand 
side  of  Eq.  (25)  to  be  evaluated  analytically. 

These  expressions  are,  of  course,  not  accurate  too  near  the  edges  of 
microstrip  patch  and  we  dare  not  let  A  =  0. 

However,  since  w/h»l,  we  can  safely  employ  near  each  edge,  a  vec¬ 


tor  potential  of  the  form  given  in  Eq.  (A-4)  .  In  particular  with, 

H  h 


(28  a) 


H'ds  =  u,  (x  =  J,  y  =  h) 


the  conjugate  potential  functions  u  and  v  satisfy 


(28  b) 


Jii-s/li  .  eii  cosv  t  u  ,  ! 


(29a) 


^  =  eu  sinv  +  v 


(29b) 


near  the  right  hand  edge  (x  *  w/2). 

The  last  integral  on  the  right  hand  side  of  Eq.  (25)  can  be  evaluated 


H£(e)  ^  (ua  •  ub> 
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where  from  Eqs.  (20a, b). 


e  a  +  u  +  1 

a 


«4  "b  , 

TT  '  e  -  ub  -  ' 


For  convenience  in  evaluating  the  first  three  integrals,  we  choose  A  =  h/n 
and  find  2(ua  -  u^)  =  -1.09.  From  symmetry  consideration. 


H  (O  d£,  =  H  (O  df. 


therefore  the  complete  analytic  approximation  of  Eq.  (25)  may  be  written  as 


M  H  (hl)e' 


Z  4ttT 


sin(ksin6cos>i>v) 


x  sin(-^  sinecos<j>)  -  cos(-k^  sin0cos4>) 


Because  kh/TT«l,  the  lower  limit  of  the  integral  can  be  set  equal  to  zero 

00 

and  since  sina£.-  a  ,  the  final  result  is 

r dt-  "/2  in 


pH  (hL)  ..  r 

a  o  o'  _-.ikr  ,-,„x  .  kw 

z - 4^ —  e  F(a)  S1n  ~  S1 


Y  sinocos*  .55  (cos  ^  3 i nOcos^  )  I 


where  h  (t>)  is  given  by  Eq.  (26). 

P  _ 

Of  course,  and  t  are  related  by  E  /H„  =  ~ — i —  with  r*  calculated 

0  0  0  0  k  eff  ro  eff 

from  Eqs.  (A-9)  and  (A-10). 

Notice  that  this  estimate  of  A?  assumes  that  the  ground  plane  extends 
many  wavelengths  and  also  that  the  transverse  field  pattern  is  Laplacian  in 
nature  (a  consequence  of  the  TEM  approximation  that  neglects  any  z-components 
of  the  field  generated  by  the  fringing  at  the  open  front  face). 

Eqs.  (16)  can  now  be  employed  to  calculate  the  radiation  fields, 
radiated  power  and  Gr«  The  latter  value  can  be  used  to  modify  Eq.  (17) 
if  necessary. 
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3.  COUPLING  BETWEEN  STACKED  MICROSTRIP  PATCHES 

In  order  to  accomodate  dual -frequency  operation.  Ball  Aerospace  Sys¬ 
tems  Division  has  developed  the  stacked  element  concept  4  in  which  a  smaller, 
high-frequency  patch  is  placed  above  the  larger,  low-frequency  patch.  In 
this  section,  we  develop  a  transmission  line  model  with  distributed  coupling 
and  solve  for  the  excitation  in  both  patches  when  the  drive  frequency  cor¬ 
responds  to  either  the  low  or  high  resonance  condition.  We  consider  the  cases 
of  substrates  both  with  dissimilar  and  with  identical  dielectric  properties. 

A  cross-section  of  the  stacked  micro-patch  geometry  is  shown  in  Fig. 
7a  and  represents  any  x-y  plane  in  the  range  0<z<L^ .  The  two  dielectric  slabs 
need  not  be  identical.  A  section  through  the  y-z  plane  is  given  in 
Fig.  7b.  A  transmission  line  model  for  the  z-variation  of  the  fields  is 
shown  in  Fig.  7c.  Notice  that  coupling  region  extends  over  the  length  of  t.he 
upper  plate,  (0  s  z  *  L). 

3.1  TRANSMISSION  LINE  MODEL 

The  transmission  line  equations  governing  the  voltages  V^z)  and  V^U), 
measured  with  respect  to  the  bottom  ground  plane,  and  the  currents  i-j(z) 
and  ipU),  representing  respectively  the  total  currents  on  the  conducting 
patches,  are  for  quasi-TEM  modes,  wel 1 -approximated  by 

3i,  3v,  3v9 

"  ~9z  =  cil  TF  +  C12  “at  (34a^ 


3V,  3i. 

I  -  I  I  I 


,  9i2 


_ L  =  I  '  L  +  I  '  £ 

3Z  L11  at  12  3t 


(34b) 
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The  high  frequency  microstrip  patch  is  placed  above  the  low 
frequency  radiator  as  shown  in  a  and  b.  The  equivalent 
transmission  line  circuit  is  shown  in  c. 

FIGURE  7.  STACKED  ELEMENTS  FOR  DUAL -FREQUENCY  COVERAGE 
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and 


(35a) 


ai2 

Cl 

9vl 

+  C '  „ 

av2 

az 

H2 

at 

c22 

at 

V  ) 

3v2 

i 

1 _ 

ail 

+  LA 

ai2 

(35b) 

IT  = 

l12 

at 

l22 

TT 

where  Cj2  <  0  and  Lj2  >  0. 

Because  these  voltages  and  currents  contain  contributions  from  the 
fields  under  both  patches  even  if  the  fringing  fields  are  zero»the  coupling 
terms  are  always  non  zero. 

It  is  therefore  preferable  to  define  new  variables,  namely 


va  =  v1  -  v2  (36a) 

vb  =  v2  (36b) 

and 

ia  =  (37a) 

ib  =  i2  +  (37b) 


The  coupled-transmission  lines  for  these  choices  are  shown  in  Fig. 
7c;  the  governing  equations  are 


9ia 

az 

-cia 

9va 

at 

+ 

c;b 

8vb 

at 

(38a) 

ava 

aia 

ai. 

a 

az 

=  L« 

at 

+ 

Lab 

b 

at 

(38b) 

aih 

ava 

avK 

b 

az 

Cab 

a 

at 

+ 

Cbb 

b 

at 

(39a) 

av. 

ai 

ai. 

b 

az 

mLib 

a 

at 

+ 

Lbb 

b 

at 

(39b) 
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where  the  capacitances  and  Inductances  per  unit  length  (along  z)  are 


1  1 

Caa  =  Cil  ' 

(40a) 

1  pW  p 

-  cj,  ♦  c-2  ♦  2C|2  - 

(40b) 

Cab  =  C11  +  C12  -  keFaaCbbs  0 

(40c) 

hl 

L1  =  L  '  +  1  •  _  ?i  •  ~  _L 

Laa  L11  l22  ^l12  *'o  W] 

(41a) 

hp 

L'  =  c- 

bb  22  o  w2 

(41b) 

i  '  =  i  1  -  i  1  =  l  i  1  i  1  >n 

Lab  l12  l22  -  KmfaaLbb 

(41c) 

The  electric  and  magnetic  coupling  coeffients,  respectively  k^  and 
are  both  less  than  unity;  if  fringing  is  neglected,  they  both  vanish  and 
the  lines  are  uncoupled. 

Decause  the  mode  Qs  are  moderate  to  high,  it  is  reasonable  to  neglect 
the  radiation  from  the  open  circuited  ends  of  the  lines,  this  implies 
that  the  mode  voltages  and  currents  are  in  time-guadrature. 

We  therefore  take 


Vb  '  Va,b(*> 

‘a.b  *  ‘a,b<z> 


and  reduce  Eqs.  (38)  and  (39)  to 


al  3  (oVa  +  “'Cab  Vb 


dVa 

dz 


Ml 

U  -  w  I  +  <oL'k  Ik 
o  w |  a  ab  b 


(42a) 

(42b) 


(43a) 

(43b) 
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(44a) 


d  I  i  t  nW  n 

^  0)  V.  +  0)C '  V 

dz  h2  b  ab  a 


dV.  h, 

~3^  =  v‘o  W~  wIb  +  wLab  la  * 


(44b) 


These  in  turn  may  be  combined  to  yield 


+  <Ai0c,  *  LibCab*  V  "“^'o  S"J  Cab  *  c2  Lab>  Vb  (45a) 

+  “2‘>‘oc2  *  Lib  cab>vb  5  -“2<"o  |  Cib  +  C1  h}  Kt>>  Va.  (45b) 

If  the  two  dielectric  constants  are  the  same,  cj  =  c2  =  e,  the  uncoupled 
phase  velocities  of  the  two  plate  transmission  lines  are  equal  and  coupling 
is  enhanced. 

In  general,  the  spatial  voltages  along  the  two  lines  are  found  from 


Eq.  (45)  to  be 


Vg  =  A  sink-jZ  +  B  sink2z  (v  s  z  s  L^)  (46a) 

k2  ‘  ka 


It  -  |< 

kl  Kb 


A  sink.z  +  B  sink?z  (0  <  z  <  L,) 


(46b) 


C  cos  kb  (L2  -  z) 


(L,  S  z  <  Lo)  . 


No  mutual  coupling  is  assumed  to  exist  over  the  interval  L-j  i  z  <  L2- 
Here,  because  the  coupling  is  weak, 

ka  =  w2(Vl  +  Lab  Cab)aa>2Vl 


(47a) 


kb  =  &)2^Moe2  +  Lab  Cab^  =a)2lloc2 


(47b) 
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o  ?  "?wl 

k  i  =  co  p  t— — —  Jr”.  (  fc  _  k  +  It  .  k  ) 

I  o  '  1 '  ’  2  e  1  1  nr 


(48a) 


and  kj  ^  are  the  roots  of 


(k2-  k2)  (k2-k2)  =  k2  k2j 


The  associated  currents  are  from  Eqs.  (43b)  and  (44b), 


J _  wlw2  /  h2  dVa  _  ,  ,  dVb^ 

■,2s  ^  h  h  \0  W9  "  ab  <*2  ' 

-  \?  >'o  hlh2  V  2 


.  _  1  "1*2  (  hl  dVb  dVa  ) 

b '  ^  w  '  4b  ■ 
m  o  1  c 


(50a) 


(50b) 


When  boundary  conditions  are  imposed  (continuity  of  Vh  and  I b  at  z  =  L ^ 
and  Ia(L^ )  =  Ib(L,>)  =  0).  the  coefficients  A,  B  and  C  are  constrained  to 


/  h2  .  ,  kII.'S\L  (\  |  +  (y  —  T  --r  -  L'A  k„  B  cos  k„l,  *0 

11  Vwc  „ 


(51a) 


kc  -  k 
K1  v 


Asink^L^  +  B  sin  k^L  •(  =  C  cos  kb(L2-  )  (51b) 


(-4b*"oij  ^jklA“sklLl 

^bTl  *uoi£)k2BcOSk2Ll  *"o»l  kbCsinkb(L2-Ll'- 
V  k2“Ka  x/  ( 


JS)k1  A  cos  k^L^ 


22 


In  order  to  establish  the  resonance  conditions,  the  determinant 
formed  from  Eq.  (51)  is  set  equal  to  zero  and  the  characteristic  frequencies 
must  be  found. 

As  expected  "x/4  -  resonances"  exist  for  either  k  L.~-n/2  or 

3  I — 

k^2  =  tt/2 .  If  the  two  frequencies  ^  and  a>2  are  too  closely  spaced,  ap¬ 
preciable  excitation  can  occur  on  both  lines.  We  next  consider  the  unwanted 
excitation  on  line  b  when  line  a  is  driven  at  its  frequency,  and  vice  versa. 
Two  situations  are  of  interest.  In  the  first  instance  Cj  /  and  the 
phase  velocities  of  the  two  lines  are  unequal;  in  the  second,  =  e  aucj 

the  velocities  are  also  equal. 

3.2  UNEQUAL  PHASE  VELOCITIES 


When  e-|  /  c2  and  the  coupling  is  weak 

+  kI 

T - IT" 


kl  =  ka 


(52a) 


and 


k2  ~  kb  * 


x2  k  2 

—L1L1 

2kb«kb  - “Ip 


(52b) 


If  line  #  2  is  driven  so  that 


VVo  si"ko2 


(59a) 


it  follows  that 


Notice  that  if  a>1  the  denominators  of  Eqs.  (53b),  (54b),  (58b) 
and  (59b)  will  all  tend  to  vanish  leading  to  large  amplitudes  of  voltage  in 
the  undriven  lines. 


Even  neglecting  loss  as  we  have  done,  the  true  amplitudes  will  not  be¬ 
come  infinite  when  the  more  exact  Eqs.  (51)  are  considered.  The  crucial  point 
however,  is  that  the  approximate  results  allow  comparisons  to  be  made  among 
the  cases  considered  and  indicate  how  close  the  dual  resonances  can  become 
before  appreciable  mode  excitation  in  the  undriven  patch  creates  problems 
with  the  radiation  characteristics  of  the  composite  antenna. 

The  critical  factor  is  Icos  n  ^2  I"1.  Notice  that  for  dual -frequency 

I  7^1 

applications  of  current  interest,  namely  the  Global  Positioning  System  (GPS) 

and  Aerosat  ,  the  values  of  are  respectively  1.3  and  1.07.  Consequently, 

the  respective  values  of  I  A  w2\  M  are  2.2  and  9.1. 

|cos\7  ay  I 
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^ •  CONCLUSIONS  AND  RECOMMENDATIONS 

The  work  described  in  Volume  I  of  this  report  allows  predictions  to  be 
made  regarding  the  non  ideal  behavior  of  rectangular  microstrip  patch  antennas 
including  the  effects  of  ground  plane  currents  and  intra-element  coupling  due 
to  fringing  electric  and  magnetic  fields. 

It  would  be  helpful,  and  the  principal  investigator  recommends,  that 
a  sufficient  number  of  detailed  calculations  be  carried  out  for  single  and 
multiple-elements  (with  a  variety  of  substrates)  to  establish  the  range  of 
validity  of  the  formulas  developed  here  by  comparing  the  results  with 
presently  available  data. 

It  is  likely  that  some  new  experiments  would  have  to  be  carried  out 
to  allow  a  full  comparison  to  be  made. 

After  such  validation,  the  design  of  antenna  arrays  with  exacting 
pattern  requirements  should  be  attempted  and  compared  with  coordinated  ex¬ 
periments. 
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APPENDIX  A 

MICROSTRIP  FORMULAS  BASED  UPON  NOVEL  APPROXIMATIONS 
A.  1  INTRODUCTION 

The  basic  microstrip  geometry  given  in  Fig.  A1  consists  of  a  dielectric 
sheet  of  thickness  h  and  dielectric  constant  e'  £  1,  clad  on  its  lower  surface 
(y=-h)  with  a  conducting  sheet  and  on  its  upper  surface  with  a  conducting 

strip  of  width  w .  Free  space  exists  above  the  dielectric  (y^O). 

14-1 7 

Although  many  workers  have  described  and  given  approximate  ex¬ 
pressions  for  the  characteristic  impedance  and  phase  velocity  of  the  quasi-TEM 
mode  that  propagates  on  microstrip  it  is  appropriate  to  re-examine  the  subject 
and  we  do  so  in  this  appendix. 

It  is  not  that  the  previous  work  is  inadequate  for  computation  of  the 
characteristic  impedance  and  wave  velocity  of  microstrip,  but  rather  that 
we  wish  to  develop  (if  possible)  a  simple  approximation  that  will  allow  cal¬ 
culations  of  field  distributions  which  in  turn  can  be  used  to  find  electric 
and/or  magnetic  coupling  between  microstrip  patches.  Fortunately,  we  are 
able  to  achieve  our  goal  and  as  a  by  product  several  new  formulas  are  de¬ 
veloped  that  should  also  prove  useful  in  microstrip  transmission  line  analysis. 
In  particular,  the  new  formulas  for  h/w  £  1  and  arbitrary  value  of  dielectric 
constant  are  particularly  simple  and  ZQ(h/w)  can  be  inverted  to  yield 
N/w(Zq),  The  latter  form  facilitates  the  synthesis  of  lines  of  prescribed 
impedance. 

A. 2  APPROXIMATE  POTENTIAL  FUNCTION 

With  reference  to  Fig.  A1 ,  and  assuming  h/w«l ,  consider  the 
function 
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(Ala) 

y  >  0, 

(Alb) 

(A2a ) 

h  s  y  <  0. 

(A2b) 

First,  it  should  be  noted  tht  this  function  was  not  derived  in  the 
usual.  Instead,  it  was  developed  as  an  interpolation  based  upon  physical 

intuition.  Second,  notice  that,  if  considered  as  a  potential  function,  Eqs. 

(1)  and  (2)  satisfy  exactly  the  boundary  condition  of  Fig.  A1 ,  namely  i^(y=-h )  =  0, 
ii,(|x|-'  w/2»  y  =  0)  =  VQ  and 

( I x |  >  w/2.  y  *  0  +)  -  e1  (|x|  >  w/2,  y  ■  0-). 

Third,  notice  that  for  [x2  +  y2  -  (|)2]2  +  (w/)2>'W2d2  and  y>0, 

y 

*l>  -  T(x,y)  which  is  a  Laplacian  function,  as  is 

♦i2  +  y2<<w  ,  y  >  o)  *  vQ(i-  ~  y) 

and 

4<(  | x | < w/2 ,  y  <0)  =  VQ(1  +  y/h). 

Although  <p(  | x | >w/2 ,  y<  0)  is  not_exactly  Laplacian,  the  discrepancies  vanish  for 
large  |x|.  Therefore,  Eq.  (Al)  is  a  comparatively  simple  formula  that  inter¬ 
polates  between  asymptotically  correct  limits. 
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The  conducting  plane  (y  =  -h)  and  strip  (y  =  0,  |x|  <  w/2) 
are  assumed  to  have  infinite  conductivity  and  zero  thick¬ 
ness;  the  slab  of  dielectric  constant  e*  >  1  and  thickness 
h  is  assumed  lossless. 

FIGURE  A- 1 .  MICROSTRIP  GEOMETRY  SHOWN  IN  CROSS  SECTION 
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The  principal  errors  are  associated  with  the  region  near  the  plate 
edges  because  this  approximation  does  not  allow  the  E-field  to  become 
infinite  at  these  locations.  This  defect  is  important  only  when  h/w  is 
not  small  compared  to  unity. 

Additional  insight  may  be  gained  by  considering  the  potential  near  the 


When  c '  =  1,  the  exact  potential  may  be  found  by  the  method  of  conformal 
transformation ;  the  well  known  result  is 


ni  =  n(y  +_A)  .  v 

h  h  v 


cot  v  +  <n 


ib'J-JlL  .  v" 

h _ 

l  si n  v 


(A4) 


where  v  =  . 

Vo 

A  comparison  between  the  approximate  and  exact  formulae  is  given  in 
Fig.  A-2 .  Notice  that  the  E  field  strength,  proportional  to  the  spacing  between 
equipotentials ,  is  given  quite  accurately  except  very  close  to  the  edge  of 
the  upper  plate  (x  =  0,  y  =  0).  Also  notice  that  the  approximate  formula 
(and  also  physical  reasoning)  indicates  that  the  extent  of  the  fringing  is 
reduced  as  e’  is  increased.  This  fact  suggests  that  the  approximate  formula 
becomes  more  accurate  in  such  cases  since  the  regions  of  non-Laplacian  inter¬ 
polation,  between  those  that  are  asymptotically  correct,  are  reduced. 
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For  example,  the  x1  intercept  of  the  approximate  equipotential  <p  =  VQ/ 2 
(for  y  =  0)  is  given  by  x'/h  =  1/tte'. 

In  the  event  that  the  conducting  striD  and  ground  plane  carry  dc 
current  and  the  dielectric  is  nonmagnetic,  an  analogous  formula  for  the  vector 
potential  can  be  given  by  replacing  ip  with  Az,  VQ  with  HQh  and  e'  by  unity. 

The  flux  per  unit  length  is  then  Flux'  =  yQAz  and  H  =  Vx(izAz).  Notice  that 
Hq  is  the  strength  of  the  (approximately)  uniform  field  within  the  parallel 
plate  region. 

A. 3  APPROXIMATE  MICROSTRIP  FORMULAS 
A. 3.1  Wide  Microstrip 

When  h/w«l ,  Eqs.  (1)  and  (2)  can  be  utilized  approximately  to  cal¬ 
culate  the  capacitance  and  inductance  per  unit  length  of  microstrip  trans¬ 
mission  lines. 

A. 3. 1.1  Capacitance  Per  Unit  Length 
For  |x|<w/2,  y  =  0+ 


*1; 


_ s  *wy  _ 

wh  +  A  [(y)  -  x2] 

C  -1 


(A5a) 


For  |  x |< w/2 ,  y  =  0- 


ip  -  Vn  (1  +  y/h) 


(A5b) 


The  total  free  charge  per  unit  length  of  z  on  the  upper  plate  is 


therefore 


w/2 


W2 


Q-  =  e 


V  e'w  dx 
o 


wh 


+  *c'[(£)2  -  x2] 


+  e 


V  h  dx 
o 


(A6a) 


-w/2 


-w/2 


and 


2e  * 

C'  -  g  *  *Z=  tanh- 

fe2  *  S 


'  t  .  w/2-  ~  A 
' 


(A6b) 
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A . 3 . 1 . 2 .  Inductance  Per  Unit  Length 


When  the  conducting  strip  carries  a  dc  current  I,  the  magnetic  vec¬ 
tor  potential  is  z-directed  with  A2/(HQh)  replacing  i|’/V0  and  setting 

»  '  =  1  in  Eqs.  (A-l)  and  (A-2).  Then,  since  H  =  Vx(TzA^),  it  follows  that 

+w/2 

1  =  I  [H  (x,y=0)  -  H  (x,y=0+)]  dx  .  (A-7a) 

-w 12  x  x 


Because  the  magnetic  flux  per  unit  length  (along  z)  is  Flux'  =  u^,  the 
inductance  per  unit  length  L'  is  given  by 


Flux' 


u 


0 


2w/ti_ 

(w)?  wh 

2  n 


tanh 


( w  )  w  ti 

2  „ 


(A-7b) 


It  is  convenient  to  define  a  dimensionless  parameter  o  =  4h/nW. 
In  terms  of o,  Egs,  ( A — 7b )  and  (A-6b)  become  respectively 


and 


Although  these  expressions  are  expected  to  be  valid  only  when 
h/w<<l ,  we  have  found,  empirically,  that  their  validity  can  be  extended  over 
the  entire  range  h/w  ^  1  by  adding  simple  correction  terms  to  the  denominators 
of  Eqs.  A-8a,  b  which  then  become 
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If  the  dispersion  induced  by  the  non-TEM  components  of  the  true 
transmission  line  fields  is  neglected,  certainly  valid  for  h/w«l ,  the 
phase  and  group  velocities  are  given  by  v  =  c/v^*'e^f  where 

e'eff  =  c  Jro7"  .  (A-10) 


The  appropriate  characteristic  impedance  is  found  from 


Zo  = 


L' 

C* 


(A— 1 1 ) 


The  accuracy  of  our  approximations  can  be  assessed  by  comparing  the  values  of 

13  15 

e'  «  and  Z  e  with  the  accurate  values  of  Wheeler  and  Hammerstad  that 
eff  o  o 

were  developed  by  a  combination  of  exact  and  quasi-empirical  means.  This  compari¬ 
son  is  given  in  Table  1  for  both  Eqs,  (A-8a,b)  and  Eqs.  (A-9a,b). 


A. 3. 2  NARROW  MICROSTRIP 

When  w/h<l,  the  approximate  formulas  given  above  are  not  valid  and 
it  is  instructive  to  consider  an  alternative  approach  based  upon  the  method 
of  images. 

First  assume  c'  =  1,  and  note  that  the  charge  distribution  of  Fig. 
(A-3)a  produces  the  proper  electric  field  in  the  half  space  y  ;>  0  provided 
Q'(x)  and  its  image  cause  the  strip  y  =  h,  |x|  S  w/2  to  be  an  equipotential 
of  value  V.  If  w/h  is  sufficiently  small,  the  effect  of  the  image  charges 
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TABLE  A-l  WIDE  MICROSTRIP 


h/w 

e'eff 

z/[5 
°  K 

Morqenthaler 

(A-8a,b)(A-9a,b) 

Wheeler 

Hammers tad 

Morqenthaler 
(A-8a ,b)(A-9a,b) 

Wheeler 

Hammers  tad 

.01 

1 

1 

.0096 

.0096 

.05 

1 

1 

.0443 

.0427 

.0427 

.10 

1 

1 

.0826 

.0767 

e '  * 

.50 

1 

1 

.3284 

.2362 

.2350 

.2348 

1.00 

1 

1 

.5972 

.3339 

.3346 

.10 

1.4033 

1.4017 

1.4050 

1 .4044 

.0069 

.0068 

.0069 

.0068 

.05 

1.3808 

1.3723 

1.3785 

1.3767 

.0321 

.0311 

.0311 

.0310 

.10 

1.3657 

1.3497 

1.3579 

1.3554 

.0605 

.0568 

.0564 

.0567 

C  '  2 

.50 

1.3419 

1.284 

1  .3030 

1.303 

.2448 

.1839 

.1804 

.1824 

1.00 

1.3485 

1.2439 

1  .2837 

1.2834 

.4429 

.2663 

.2610 

.2607 

.01 

2.2067 

2.2039 

2.2111 

2.2113 

.0044 

.0044 

.0044 

.0043 

.05 

2.1427 

2.1286 

2.1385 

2.1404 

.0207 

.0201 

.0200 

.0199 

.10 

2.0963 

2.071 

2.0816 

2.0853 

.0394 

.0370 

.0368 

.0368 

c'  = 

.50 

1.9938 

1.9183 

1.9308 

1.9380 

.1647 

.1231 

.1217 

.1223 

1.00 

1.9901 

1.8718 

1  .8799 

1.8854 

.3001 

.1784 

.1782 

.1775 

.01 

3.1142 

3.1110 

3.1215 

3.1228 

.0031 

.0031 

.0031 

.0031 

.05 

3.0068 

2.9920 

3.0024 

3.0096 

.0147 

.0143 

.0143 

.0142 

.10 

2.9262 

2.9019 

2.9080 

2.9213 

.0282 

.0264 

.0263 

.0263 

r'  «  ' 

.50 

2.7230 

2.6795 

2.6612 

2.6834 

.1206 

.0881 

.0883 

.0884 

1.00 

2.6886 

2.6348 

2.5792 

2.5977 

.2221 

.1267 

.1299 

.1288 
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on  the  form  of  Q'(x)  is  negligible  and  the  latter  can  be  found  from  the 
solution  of  a  single  charged  strip.  That  potential  is  proportional  to  v 
(taken  as  zero  on  the  strip)  where 

- +  =  (?)2  (A_12) 

1+sinh  v  sinh  v 


The  potential  that  is  the  superposition  of  a  pair  of  such  terms  (the  second 
due  to  the  image)  is 


♦  *  iSr  sinh‘1 
0 


i  *  «<yj->2  - 

1 

- 

+  \[(2~)2  +  4(^-h)2  -  l] 

n|  2 

r(~)2  +  4(^)2-l' 

+ 

\[£)2  *  Hylh)2  -  1 

2 


(A- 13) 


Although  <H|x|  /2,  y  =  h)  should  be  exactly  constant,  the  approxima¬ 

tion  is  not  precise  but  for  simplicity  set  <Ho,h)  =  V.  Then 

V  =  sinh"1  (^h)  .  (A- 14) 

0 

Since  Q*  =  C'V,  it  follows  that 


0  1  ~  1  sinh’1  (~) 

O'  "0  271  w 


(A-  IS) 
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4 _ -  a  yoii 


IMBH 


When  f.'  >1,  the  situation  is  more  complicated,  but  the  method  of  images 
can  still  be  used. 


The  image  pair  with  dielectric  slab  of  thickness  2h  between  them  showr 
in  Fig.A3a  has  an  E-field  within  the  dielectric  that  is  exactly  equivalent  to 
the  free  space  problem  with  an  infinite  set  of  images.  The  kth  set  of  images 
is  given  by  +  q^x)_2___  ^e'-l^  k  Images  for  k  =  0,  1,  and  2  are  in¬ 

dicated  in  Fig.  A-3b. 


If  h/w  is  large  enough  so  that  the  noninteracting  strip  approximation 
can  be  rnade.it  follows  from  superposition  and  evaluation  of  g/( x=0 ,  y=h)  =  V, 
that 


Lo  1 

r-TP+niF 


sinh'1^)  +  T.  (fe‘) 


•  ir  r 


w 


k=l 


1+c ' 


si 


Jnh"1(--(k+lj)  -  sinlT-1(^k) 


Because  by  assumption  £  >  1 ,  the  values  of 


sinh'1  (~  n)  =  nn 


£n*/(£n)2  +  l 


where  n 
Then 

and 


k  or  k  +  1  ,can  be  reasonably  approximated  as  K,n(~n  ) 
sinh-1  Q  (k+1 )  -  sinh-1  k)  -  n(l  +1) 

jsinh-1  {—)  ♦  (]^-)k  tn(l+J)j 


(A-16) 


(A-17) 


( A- 18) 


(A- 19) 


Notice  that  the  geometric  factor  h/W  has  disappeared  from  the  infinite  sum. 
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)'  --  ('  \  o  ' 

•fk  .  -.1  \r^l)  % 


Fields  in  the  upper  half  of  a  symmetric  structure  consisting 
of  a  double  thickness  dielectric  slab  and  a  second  conducting 
strip  carrying  image  electric  charges,  as  shown  in  a,  are 
equivalent  to  those  in  Fig.  A-l.  Insofar  as  the  E-fieid 
inside  the  dielectric  is  concerned,  the  dielectric  can  be 
replaced  by  the  image  of  density  Qj^  shown  in  b. 


FIGURE  A-3.  REPLACEMENT  OF  MICROSTRIP  DIELECTRIC  LOADING 
WITH  IMAGE  CHARGES 


L ^ ^ ^  i.  __  w- 


For  convenience,  we  define  a=  -i-^y  and  since  a  lies  between  zero  and  unity 
it  is  possible  to  approximate  the  summation  with  a  low-order  polynomial.  A  very 
simple  and  quite  accurate  approximation,  is 


oo 

I  (-a)  tn(l+£-)=  -  tHa-®  (23- 12a  +  4a2)  =6 (a) 


k:1  '  '  ' '  k  TOO  1 


valid  for  all  0<  a  <  1 

The  final  set  of  equations  valid  for  h  >  1  is 


These  are  simple  and  lend  themselves  to  inversion  since 


\/¥z°  Vo  (?) 


^  Zo  ♦  <W2>2"  -  6/2)  2 

ft  -  X  (e'-l)  / 23-12(e'-l )  .  we'-l)2) 

6  m  r+r  v  prr  +  4(^+t  / 


(A-20) 


L’  =  1  sinh'1  (4h) 

M0  2tt  w 

(A-21a) 

eo  _  1 

Tsinh"1  —  +  3  ( e '  - 1 )  1 

C  (e  +1  )ir 

L  w  e'+l  J 

(A-2 1 b ) 

(A-22) 


(A-23) 


(A-24) 


Neglecting  dispersion,  the  slowing  factor  c/v  =  e'  is  again  given  Dy 
Eq.  (A-10) .  Numerical  comparison  with  the  published  approximate  formulas 
of  Wheeler14  and  Hammerstad16  is  very  favorable,  as  shown  in  Table  A-2. 


TABLE  A-2  NARROW  MICROSTRIP 


& 

eff 

vf* 

h/w 

Morgenthaler 

Wheeler  Hammers tad 

Morgenthaler 

Wheeler  Hammers tad 

Eqs.  A-21a,b 

Eqs.  A-21a,b 

1 

1 

1 

.3334 

.3350  .3359 

2 

1 

1 

.4419 

.4423  .4425 

5 

1 

1 

.5872 

.5873  .5873  c'  *  1 

10 

1 

1 

.6974 

.6975  .6975 

100 

1 

1 

1 .0639 

1.0639  1.0639 

1 

1.2837  1.2801 

.2593 

.2610  .2624 

2 

1.2700 

1.2701  1.2669 

.3479 

.3483  .3493 

5 

1.2584 

1.2589  1.2557 

.4666 

.4665  .4677  e'  =  2 

10 

1.2529 

1.2534  1.2497 

.5567 

.5564  .5581 

100 

1.2430 

1.2434  1.2386 

.8859 

.8556  .8590 

1 

1.8799  1.8854 

.1764 

.1782  .1781 

2 

1.8473 

1.8450  1.8493 

.2392 

.2397  .2393 

5 

1.8166 

1.8165  1.8186 

.3232 

.3233  .3229  c'  =  5 

10 

1.8024 

1.8027  1.8018 

.3869 

.3869  .3871 

100 

1.7778  1.7709 

.5986 

.5984  .6008 

1 

2.5991 

2.5792  2.5977 

.1283 

.1299  .1293 

2 

2.5292 

2.5234  2.5387 

.1747 

.1753  .1743 

5 

2.4797 

2.4781  2.4882 

.2368 

.2370  .2360  e'  »  10 

10 

2.4569 

2.4563  2.4607 

.2839 

.2840  .2834 

100 

2.4166 

2.4169  2.4096 

.4402 

.4402  .4415 
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A. 4  QUASI-STATIC  COUPLING  BETWEEN  MICROSTRIP  PATCHES 

As  shown  in  Fig.  A4,  consider  the  two  microstrip  patches  of  widths 
Wj  and  w2  separated  by  a  distance  l.  The  substrate  has  a  dielectric  constant 
e'  and  thickness  h. 

It  is  well  known  that  the  transverse  field  patterns  associated  with 
a  TEM  mode  are  Laplacian  in  character  and  are  therefore  identical  in  form 
to  those  generated  by  dc  voltages  or  dc  currents.  It  is  this  fact  that 
allows  the  use  of  quasi -statics  to  calculate  electric-and  magnetic-field 
coupling  coefficients. 


A. 4.1  ELECTRIC  FIELD  COUPLING 

If  the  patches  are  maintained  at  electrostatic  potentials  and  V2  with 
respect  to  the  ground  plane  potential  <p  =  0,  the  charges  (per  unit  length  along 

z)  on  the  two  patches  are  given  by 

Qj  =  CJ,  V,  *  c;2  V2  (A-25a) 


Q2  =  C21  V1  +  C22  V2 


(A-25b) 


where  the  mutual  capacitances  satisfy  C12  =  C21  <  0  and  it  is  usual  to  define 


a  dimensionless  coupling  coefficient  as 


Cl 2  C21 
£Tl  C22 


(A-2b) 


If  V2  =  0  and  h«  t,  it  follows  that  the  potential  below  patch  2  is  ap¬ 
proximately  zero  and  that  above  aoproximately  given  by  Eq.  (Al-a). 


Since 


l  +  w2  +  Wj/2 


Ey(x,  y  =  0+)  dx 
i  +  w^/2 


(A-2 1  a) 
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'>1 


-h 


■w,/. 


V: 


««V’ 


{♦Wj  *w, 


Two  microstrip  patches  of  width  Wj  and  w7 ,  both 
spaced  a  distance  h  above  the  same  ground  plane, 
parallel  to  one  another  and  separated  by  a  distance  2.. 

FIGURE  A- 4 .  COUPLED  MICROSTRIP  PATCHES  OF  FINITE  WIDTH 
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and 


Ey(x,  y  =  0+)  dx 


-w,/2 


the  ratio  yields 


°2  I  °2J 
"cn- 


U+  y  +w2  '  pl)  U  +  y  +p,) 
(e+  1  +w2  +  P])  +  -  p.) 


^  +  tanh-1  (!lZ2) 

h  7Tq1  q]  ' 


On  the  other  hand,  if  =  0, 


i-Pi> 


(A-27b) 


(A-28) 


(A-29) 


gives  a  syrmnetric  formula  in  terms  of  an  interchange  of  w1  and  w2  and  p2,q2  where 
I  wi  ?  wih  I  w.  „  wTh 

Pi  =  '  iP"  ’  qi  =  4(T}  +  Hi* 


w.  2  wi h 

Since  (y)  »— .  .  an  approximate  formula  for  |k|,  the  coupling  coefficient,  is 


where  L]2  -  and  k2  =  f1 


An  analoqous  analysis  with  A  instead  of  w 
1T22  2 

and  uQAz  =  gives  the  same  formula  for  |k|  as  Eq.  (A-30)  except  that  e'  is 
replaced  by  unity. 

The  magnetostatic  fields  near  the  coupling  region  when  w^  2»£  can 
be  described  by  allowing  +  “  and  using  the  method  of  conformal  trans¬ 
formation. 

For  the  boundary  conditions  shown  in  Fig.  (A5),  the  H-field  can  be  ob¬ 
tained  fromvx  A  where  Az  =  Cv,  H  •  ds  =  Cu  and  the  conjugate  potential 
functions  are  related  by 


ir(b2  -  l)x  _u  ,  ,  b2  -  1 

=  e  cosv  -  1  + - ^ — 


tn 


4  -  4e  cosv 


(A- 32a) 


and 


rr(b  -  1  )y 
2h 


e  sinv 


tan 


Z'  2sinv  N 
\2cosv  -  ey 


(A-32  b) 


The  constant  b  is  the  solution  of  the  transcendental  equation 

2b  „„/b+l '  _  ttS 

^  +  W(PT>  -  2h 


(A- 33) 


A.4.3  COUPLING  BETWEEN  A/4  RESONANT  PATCHES 

Reference  to  Fig.  A6a  reveals  that  the  mode  patterns  dictate  that 

magnetic  coupling  will  dominate  near  the  shorted  end,  whereas  electric  coupling 

will  dominate  near  the  open  circuited  end. 

Patches  oriented  as  shown  in  Fig.  A6b  have  very  weak  H-fields  near  the 

« 

open  circuited  edges,  hence,  only  electric  field  coupling  is  important. 

In  all  cases,  an  increased  dielectric  constant  reduces  the  fringing  E 


fields,  hence,  the  value  of  i^eiectr^cl  by  a  factor  that  approaches  e'. 
(Note  that  In  this  case 

the  coupling  formulas.) 


(Note  that  In  this  case  and  L2  should  be  used  instead  of  w-j  and  w2  in 
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Identical  to  Fig.  A-4,  except  that  =  w9  ■><».  The 
static  electric  and  magnetic  fields  generated  by  dc 
charges  and  currents  are  solved  by  the  method  of 
conformal  mapping. 


FIGURE  A- 5.  COUPLED  MICROSTRIP  PATCHES  OF  INFINTE  EXTENT 
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Two  rectangular  patch  microstrip  antennas  placed  either 
side  by  side  (displaced  a  distance  H  along  x)  or  in  line 
(displaced  a  distance  l  along  z)  both  electric  fields 
and  magnetic  fields  couple  the  patches  shown  in  a),  whereas 
it  is  principally  electric  fields  that  couple  b) .  Notice 
the  location  of  the  short  circuited  edges  in  b. 

FIGURE  A- 6 .  COUPLED  RECTANGULAR  PATCH  MICROSTRIP  ANTENNAS 
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A* 


Once  | k |  is  known,  the  amplitude  of  the  mode  induced  in  patch  #  2  by 
a  known  mode  amolitude  in  patch  #  1  can  be  estimated. 

For  example,  for  the  case  of  Fig.  (A6b),  Q£  =  0,  both  before  and  after 
the  voltage  is  applied.  Therefore,  from  (A-26) 


'C21V1  +  C22V2  =  °”  ”  (A- 34a) 

Qj  -  C^(l  -  k2)V1  (A - 34 b ) 


and 


-C‘ 

V2(z)  =  V ( z )  =  |k 
c  l22  M 


41  v  (Z)  .|k 

l22  1 


J5- 


(2) 


(A- 35) 


The  coupling  between  patches  implies  that  near  resonance,  an  equivalent 
circuit  can  be  made,  of  the  type  shown  in  Fig.  A7  where  mutual  inductance  M 
is  present,  in  addition  to  the  coupling  capacitance  Cc-  The  transformers  model 
the  coupling  (via  microstrip)  into  the  resonant  natches. 

If,  as  is  common,  the  two  patches  are  identical,  the  coupling  splits 
the  degeneracy  into  closely  spaced  frequencies. 

For  circularly  polarized  arrays.  Ball  Aerospace  Systems  Division  has 
developed  a  crossed-slot  element  consisting  of  four  X/4  microstrip  patches 
arranged  as  shown  in  Fig.  A-8.  The  type  of  analysis  outlined  above  can  also 
be  helpful  in  estimating  the  coupling  coefficients.  When  known,  the  level  of 
rf  excitation  induced  on  patch  #  2  by  the  level  on  #  1  (and  vice  versa),  can 
be  estimated.  Then  the  complete  A  and  A*  can  be  computed  which  in  turn  yield 
the  far-field  radiation  pattern  that  includes  the  effects  of  the  coupling. 
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Equivalent  circuit  governing  the  coupled  resonators  of 
Fig.  A-6.  The  mutal  inductance  M  accounts  for  the  H-field 
coupling;  the  capacitance  Cq  is  due  to  E-field  coupling. 


FIGURE  A- 7 .  EQUIVALENT  CIRCUIT  GOVERNING  COUPLED  RESONATORS 
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short  circuit 


Crossed-slot  element  consisting  of  four  rectangular 
microstrip  patches.  When  excited  properly,  circularly 
polarized  radiation  fields  are  created. 


FIGURE  A-8 .  CROSSED- SLOT  ELEMENT  FOR  PRODUCING  CIRCULAR- 
POLARIZED  RADIATION 
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APPENDIX  B  "  REPORT  OF  NEW  TECHNOLOGY 


There  are  no  inventions  or  patentable  items  contained  in  Volume  I  of 
this  report.  However,  certain  advances  in  the  theory  of  microstrip  trans¬ 
mission  lines  and  antennas  are  reported  here. 

In  particular,  for  rectangular  patch  antennas,  we  include  the  effects 
of  currents  induced  on  the  ground  planes  by  fringing  fields. 

Coupling  between  two  patches  sharing  the  same  substrate  and  ground 
plane  or  else  employing  separate  ones,  stacked  vertically,  is  also  considered 
by  means  of  a  novel  approximation  that  helps  provide  physical  insight  with 
respect  to  field  patterns,  coupling  between  patches  and  the  like. 

As  an  important  by-product  of  this  work,  several  new  approximate  formu¬ 
las  are  obtained  that  very  accurately  predict  the  electrical  characteristics 
of  microstrip  transmission  lines  of  arbitrary  width  and  substrate  thickness 
when  the  dielectric  constant  of  the  substrate  is  also  arbitrary. 

Volume  II  of  this  report  is  devoted  to  the  analysis  and  synthesis  of 
multi-resonant  elements  with  emphasis  on  dual -frequency  operation  of  rec¬ 
tangular  microstrip  patch  antennas  with  or  without  external  matching  networks 
and  mav  contain  inventions  or  patentable  items.  For  a  discussion,  please  refer 
to  the  appendix  of  Vol.  II  of  this  report. 
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